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The spectral difference between solar flare HXR coronal and 
footpoint sources due to wave-particle interactions. 
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ABSTRACT 

Aims. Investigate the spatial and spectral evolution of hard X-ray (HXR) emission from flare accelerated electron beams subject to 
collisional transport and wave-particle interactions in the solar atmosphere. 

Methods. We numerically follow the propagation of a power-law of accelerated electrons in ID space and time with the response of 
the background plasma in the form of Langmuir waves using the quasilinear approximation. 

Results. We find that the addition of wave-particle interactions to collisional transport for a transient initially injected electron beam 
flattens the spectrum of the footpoint source. The coronal source is unchanged and so the difference in the spectral indices between 
the coronal and footpoint sources is Ay > 2, which is larger than expected from purely collisional transport. A steady-state beam 
shows little difference between the two cases, as has been previously found, as a transiently injected electron beam is required to 
produce significant wave growth, especially at higher velocities. With this transiently injected beam the wave-particle interactions 
dominate in the corona whereas the collisional losses dominate in the chromosphere. The shape of the spectrum is different with 
increasing electron beam density in the wave-particle interaction case whereas with purely collisional transport only the normalisation 
is changed. We also find that the starting height of the source electron beam above the photosphere affects the spectral index of 
the footpoint when Langmuir wave growth is included. This may account for the differing spectral indices found between double 
footpoints if asymmetrical injection has occurred in the flaring loop. 
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1. Introduction 

Solar flares are prolific accelerators of electrons, sending beams 
outwards into interplanetary space and downwards deeper into 
the solar atmosphere. The latter is inferred from radio and hard 
X-ray (normally above > 20 keV, HXR) observations which 
suggest coronally accelerated electrons propagating downwards 
along magnetic field lines towards the denser chromosphere, 
eventually s topping via Coulomb colhsions heating the local 
plasma (e.g. ' Brownl 19711: iBrown et~al ]|200?). This heated chro- 
mospheric material expands upwards filling the magnetic loops, 
producing bright emission in soft X-rays SXR and EUV. In 
this standard model a power-law of accelerated electrons in the 
corona Fo(E) oc E^''*' only encounters Coulomb collisions with 
the background plasma whilst it propagates to the chromosphere. 
Assuming a time and spatially independent electron injection, 
the footpoint source will emit HXR with spectrum /(e) oc e'^p 
where ypp - - 1, namely the t hick target model (iBrownl 
[19711 ISvrovatskii & ShmelevalfTOTl . The same beam of accel- 
erated electrons should continuously emit via bremsstrahlung 
with the background plasma as it propagates from the corona. 
This t hin- target emission (assuming no energy loss by the elec- 
trons) ('Amol dv et al]ll968t iHoIt & Clinelll968t lTakakurallI969t 
[Un & Hudsonl ll976h is again found to produce a power-law 
spectrum but this time with an index of yes = 5b + 1- So if 
the HXR emission of the electron beam was observable from the 
corona to the chromosphere then the acceleration and transport 
processes involved could be well constrained. 



However the footpoint emission dominates observations due 
to the properties of the bremsstrahlung emission process, namely 
an electron distribution F{E, r) produces a photon flux emission 
/(e), in units of photons cm"^ keV"' s"', of 



/(e) cx 



(r)F(E, r)Q(e, E)dEd\ 



(1) 



where «(r) is the background plasma density and Q(e, E) is th e 
bremsstrahlung cross-section (iKoch & Mot 3[l95llHaui[l993). 
Given that the chromospheric density is several orders of mag- 
nitude greater than the corona this emission dominates the ob- 
served spatially integrated HXR spectrum. HXR imaging spec - 
troscopy i nstruments, su ch as Yohkoh/HXT (iKosugi et alJl992) . 
RHESSI ( iLin et aL '■2002), have a limited dynamic range and so 
the faint coronal emission is rarely observed with the bright foot- 
points. However observing both would be very effective for un- 
derstanding the physics of transport and acceleration of electrons 
in flares. 

The standard flare model predicts the difference between the 
observed spectral indices of the coronal and footpoint sources 
should be 



Ay = yes - ypp = 2. 



(2) 



Send offprint requests to: Hannah e-mail: iainOastro . gla .ac.uk 



While there are solar flare events with the predicted spectral 
index difference the majority of observations do not demon- 
strate this behaviour The first HXR observations of coronal 
and footpoint sources were made by Masuda et al. (1994) with 
Yohkoh/HXT and the "coronal source above two footpoints" 
geometry of this event has dominated subsequent models. The 
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spectral indices of the coronal and footpoint sources were found 
by taking the ratio of images in 14 - 23 keV to 23 - 33 keV and 
23 - 33 keV to 33 - 53 keV ene rgy bands finding y pp - 2.0, 4.0 
and yes - 2.6, 4.1 respectively dMasuda et al.lll995i) . Thus find- 
ing virtually no difference between the coronal and footpoints 
sources and certainly not Ay = 2. This event was further anal- 
ysed and found to be inconsistent with thin-target coronal emis- 
sion and thick-target footpoints (Alexan der &_Metcalf 1997J . 
Using an improved analysis technique iMasuda et al.l (l2000l) 
found that the difference in spectral indices was larger than pre- 
viously thought but still Ay < 2. They also suggest such coronal 
emission may not be non-thermal in origin but possibly due to a 
very hot (lOOMK) thermal source. RHESSI observations of such 
"above-the-loop-top sources" demonstra te power-law and not 
thermal spectra (e.g. ' Krucker et al.ll2010l) . A survey of Yohkoh 
flares with simultaneous coronal and footpoint emission found 
that the coronal emission was typically softer yet the a verage 
difference was only Ay = 1.3 + 1.5 (Petrosian et al.ll2002l) . 

With RHESSI, spatially resolved spectroscopy with high en- 
ergy resolution became routinely available and triggered a num- 
ber of interesting findings (Emslie et al. 2003; Battaglia & Benz 
[2007; Saint-Hilaire et al. 2008; B attagha & Ben z 2008; Suet all 
1200 9). Battaglia & Benz (2006) analysed five flares obtaining 
HXR spectra of simultaneous footpoint and coronal sources find- 
ing a range of < Ay < 6. Further investigation of three of 
these flares with clearly separated footpoint and coronal sources 
found that Ay > 3. Their explanation for the discrepancy is 
that there was an anomalously large density concentrated at 
the coronal source so that this emission is a combination of 
thin and thick-target, which would result in a flatt er footpoint 
spect ra l index at low energie s and A y > 2 (Bat taglia & Bend 
l2007h . IWheatland & JMelrosd (1 19951) had previously suggested 
this solution to the coronal sources observed by Yohkoh/HXT 
dMasuda et al . 1994; Feldman et al. 1994). However this model 
predicts a break in the coronal spectrum between the thin and 
thick target range which has not been observed with RHESSI 
despite its high energy resolution. A further suggestion was the 
inclusion of non-collisional losses implemented via a deceler- 
ating electric potential (Battaglia & Benz 2008). This electric 
field was assumed to be due to the return current generated by 
the electron beam and was implemented in a time-independent 
manner. It was found that this did flatten the footpoint spectrum, 
providing the necessary larger difference between spectral in- 
dices, however the simple empirical implementation had prob- 
lems with the relative brightness of the sources, in one event 
producing considerably fainter footpoint sources than observed. 

An additional constraint can be gained by considering the 
difference in spectral indices between the footpoint sources 
themselves (normally two footpoints are observed). In the stan- 
dard flare model two similar footpoints are expected to have sim- 
ilar spectral indices. However observatio ns with R HESSI have 
shown this not to be always the case. Ems lie et al.l (2003) anal- 
ysed a large X-class flare and found that during the event the 
two time correlated foo tpoint sources differed in spectral index 
with Aypp ~ 0.3 - 0.4. iBattaglia & Bend (l2006l) found that the 
difference between the footpoints' spectral index ranged from 
< Aypp < 1. Saint-Hilaire et al. (2008) analysed 53 flares 
with a double footpoint structure and found that the spectral in- 
dex between them was < Aypp < 0.6. An explanation for 
these results is that collisional transport encounter ed different 
column densities in each leg of the flaring loop (Emsl ie et al.l 
l2003h . however[siint-Hilaire et al. (2008) demonsti-ated that tiiis 
would result in the softer footpoint being brighter, the oppo- 
site to observations. They instead prefer a magnetic mirroring 
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Fig. 1. The background plasma density n(x) as a function of 
height above the photosphere. The vertical lines indicate the dif- 
ferent spatial region the emission is summed over to produce the 
total, footpoint and coronal sources. The short grey vertical line 
X() indicates a starting centroid height of the electron distribution 
at 40Mm. 



model dMelrose & Brownl Il976l) in which the magnetic field 
converges more quickly at one of the footpoints. This conver- 
gence of the magnetic field lines has been foun d to be consis- 
tent w ith area measurements of HXR footpoints dSchmahl et al.l 
I2OO6I) and has been observed in one M-class limb event with 
RHESSI llKontaret aUlIOO^ . 

In both cases the differences in spectral indices is not what is 
expected from purely collisional transport and something more 
is needed. In this paper we consider the generation of Langmuir 
waves via wave-particle interactions as the response of the back- 
ground plasma to the propagating electron beam, in addition 
to Coulomb collisions. The observation of downward moving 
decimetric radi o emission in some flares - Reverse-Slope Type 
III bursts - (e.g. T^nstrom & Zehntnerll975i Aschwanden et all 



ll995l;lKlein et alj |1997; Aschwanden & Benzlll997h is fliought 
to be a signature of Langmuir waves being generated by the 
downward moving electron beam. The presence of Langmuir 
waves is inferred as they are the starting point for non-linear pro- 
cesses that result in this non-thermal radio emission. Al though 
the majority of flares do not have such emission (Ben z et al.l 
2005) it does not mean Langmuir waves are not present but that 
the radio emission itself has been absorbed . Previous studies of 
this fast non-collisional process in flares have mostly focused 
on time and spatially independent solutions (Emslie & SmithI 
1984; Hamilton & Petrosian 1987; McClements 1987). These 
steady-state studies found that while the distribution of electrons 
changes due to the generation of Langmuir waves, the spatially 
integrated spectra shows insignificant changes over purely colli- 
sional transport. 

However, a highly intermittent and bu rsty injection of elec- 
trons is more realistic and often observed ( Kiplinger et al.ll 19841 : 
[Fleishman et all 1 19941: lAschwanden et al.lll998h and spatial fil- 
amentation (non-monolithic multi-thread flaring loops) can ex- 
plain simultaneously the heights and the sizes of observed HXR 
footpoints (Kontaretal. 2010). So to be able to simulate the 
propagation of such a possibly fragmented (temporally and spa- 
tially) electron acceleration it is crucial to consider the time de- 
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Fig. 2. Mean electron flux spectrum (left) and X-ray spectrum (right) of continually injected electron beam simulation after reach- 
ing equilibrium (steady-state), for the footpoint (orange, green), coronal (black) and total source (blue, purple). Shown are the 
simulations considering only collisional transport of an electron beam (Beam Only, dashed lines) and the addition of wave-particle 
interactions (Beam &■ Waves, solid lines). The vertical dotted lines indicated the energy range a power law was fitted over. The fitted 
spectral indices are shown in brackets (Beam Only, Beam & Waves). 



pende nt injection and transport of electrons. As iHannah et al.l 
(l2009l) have previously shown in this situation the generation 
of Langmuir waves can have a strong aff'ect on the electron dis- 
tribution and the observed HXR spectrum. 

In ^we describe the system of equations we are numeri- 
cally solving to follow the propagation of the accelerated elec- 
trons and the response of the background plasma in the form 
of Langmuir waves. In ^3.11 we demonstrate that the inclusion 
of wave-particle interactions of a continually injected electron 
beam that has reached a steady-state has little effect on the elec- 
t ron d istribution, as has been previously noted by other authors 
([Hamilton & PetrosianI 119871) . In |32] and ESwe show that 
for a transiently injected beam the wave-particle interactions do 
have an effect and show how this changes the difference in spec- 
tral indices between the coronal and footpoint sources. In ^3.51 
we show that the starting height of the beam effects the growth 
of Langmuir waves and the hardness of the resulting footpoint 
spectrum. 



2. Electron Beam Simulation 

We describe the ID electron velocity (v » vy » v^) trans- 
port of an electron beam f(v, x, t) [electrons cm"'* s] from 
the corona to chromosphere, self-consistently driving Langmuir 
waves (of spectral energy density W (v, x, t) [erg cm~^1) using the 
equations of quasi-linear relaxation (Vedenov & Velikhov" 19631 
Drummond & Pines 1964; Ryutov 1969; Hamilton & Petrosia^ 
1 9871: lKontaill2001 at IHannah et alJ l2009). This weakly turbulent 



description does not require resolving smallest scales compara- 
ble to the plasma period/Deb ye length (c.f. PIC simulations, e.g. 
iKarUckv & Kasparoval 12009 ). so the electron distribution can 
be followed at a variety of realistic heights in the solar atmo- 
sphere. We consider several processes in addition to the quasi- 
linear wave-particle interaction - Coulomb collisions. Landau 



damping and spontaneous emission - resulting in our quasilin- 
ear equations becoming 



dt dx 



Ajf-e^ d 
m^ dv 



W_df_ 

V dv 



dW 3vl dW 



dt 



dx 



n dv 



+ rcf 



d_lf 

dv 



2rL W + 



(3) 



(4) 



where n(x) the background plasma density and a)p(x)^ = 
47Tn(x)e^/m is the local plasma frequency. The resonant inter- 
action between the electrons and Langmuir waves, aip^ - kv, 
is handled by the quasilinear term, the first term on the right- 
hand side of Equations [3] and |4] The Coul omb collisions terms 
for the electrons is ycp - 47re^nlnA/m^ (iLifshitz & Pitaevskiil 
II98U lEmshd Il978h and waves jc ^ I / 3j/2Mm__^ 



Tie n]nA/(m'^Vj) 



3^ Cifshitz & Pitaevskiil [19811 iMekosd WSdti . 
In (8 X lO^n^'^^r) is the Coulomb logarithm, T 



Where In A 

is the temperature of the background plasma (taken as IMK) 
and vt = V^bTT'w is the velocity of a thermal electron, 
is the Boltzmann constant. Equation [3] without the quasilinear 
term describes the propagation of the electrons subject only 
to Coulomb collisions with the background thermal plasma. 
In Equation [4| there is also the Landau damping rate yL = 
Vi/Swp (v/vt)^ exp(-v2/2v2) jLifshitz & Pitaevskiill98lh . The 
spontaneous emission rat e is given by 5 = tj^mv In (v/vt) I (4nn) 
which agrees with that of MeIrose| ([l980l) : Pfsvtovich & Terhaan 
(ll995l) : lHanulton & PetrosianI (1 19871)" 

The evolution of energetic electrons is considered in the en- 
ergy domain above 3 keV, which is the low energy limit of 
RHESSI. The initial electron distribution is taken to be a bro- 
ken power-law in velocity which is flat below the break vq = 
3.38 X 10''cms"' (£'c(vc) = 7keV) and above it with index of 
26b (hence spectral index of 6b in energy space) i.e. 



ft , I ^ 1 / 1 if V < vc 

f(v,x,t^ 0) oc «bexp|-- ^ 



(5) 
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Fig. 3. Mean electron flux spectrum (left) and X-ray spectrum (right) of a transient initially injected electron beam simulation over 
1 second, for the footpoint (orange, green), coronal (black) and total source (blue, purple). Shown are the simulations considering 
only coUisional transport of an electron beam (Beam Only, dashed lines) and the addition of wave-particle interactions (Beam & 
Waves, solid lines). The vertical dotted lines indicated the energy range a power law was fitted over. The fitted spectral indices are 
shown in brackets (Beam Only, Beam &■ Waves.) 



where n\, is the beam density, d is the characteristic spatial 
Gaussian size. The simulation grid in x-space starts at 52Mm 
down to 0.3 Mm above the photosphere. Over this range the 
background density is at a constant 10'°cm"^ in the corona 
which then steeply rises at the transition region, exponentially 
increasing through the chromosphere, see Figure [1] The elec- 
tron distribution is initially centred at a height of ji:o=40Mm with 
d - 2 Mm. In the simulations beam densities of n^ - lO^cm"-' 
and «b - 10**cm"-' are used, meaning that the initial number 
of electrons in this simulation is ^ 23n\,d^ = 10^'"' - lO"'"* 
electrons. This is an approximation as we only have one spa- 
tial dimension to estimate the volume from, taking this as the 
FWHM - 2 \l\n2d. Compared to observations this is a mod- 
est number of electrons to be accelerated, typical for a B or C- 
class flare, however multiple bursty injections of such beams 
could achieve the numbers derived for the largest events. The 
initial wave spectral energy density is calculated from the ther- 
mal background level. 

In velocity space the simulation grid goes from Vmin - 7v't 
up to a maximum of I.OSvq, vq the maximum beam velocity. 
For the simulations shown vq = 75vt and so all over veloci- 
ties are v < c. Previous work on electron beam transport is pre- 
domina ntly non-r elativistic de spite i ncluding velocities greater 
than c ( Hamilton & Petrosianl ll987h or energies up to infinity 
dSaint-Hilaire et al.i.200 9). Wifli our upper velocity limit we re- 
main in the non-relativistic regime but we will return to this later 
in|33J] 

We use a fini te difference method to numerically solve 
Equations[3]andl4l (lKontail2001bt) . This code is written in a mod- 
ular way so that we can turn off transport processes and easily 
consider either electron transport subject only to Coulomb colli- 
sions or with the additional wave-particle interactions. For col- 
lisional transport only we solve on Equation |3] but without the 
first term on the righthand side. As there are no waves this is the 
Beam Only simulations. Considering both the collisional trans- 
port and wave-particle interactions we solve both Equations [3] 
and ID and this is termed the Beam and Waves case. This mod- 
ular nature also means that we can use an adaptive time-step in 



our code depending on the processes being used. The quasilin- 
ear relaxation is predominantly the fastest process, occurring on 
a time-scale of tq = n/in^ajp) ^ 2x 10"^ y/n/ni, seconds, where 
as the Coulomb collisions is tc = l/yc ~ 1.5 x sec- 
onds. For beam density n^ = lO^cm"-' the background density 
would need to be n > 10''*cm"^ for the collisional time-scale to 
be bigger than the quasilinear relaxation, a region electrons with 
sub-relativistic energies rarely reach. 

With these simulations we present two schemes of injecting 
the initial electron distribution. First, shown in 33.11 we have the 
continuous injection of the initial electron beam, this acting as a 
boundary condition until we reach a steady-state after about 1.5 
seconds simulation time. The other is a transient initial injection 
of the electron beam which is not replenished, shown in 33.21 
which we numerically follow for 1 seconds in simulation time, 
so that it will lose energy, completely leaving the simulation grid. 

The electron beam will be instantaneously emitting X-rays 
1(e) via bremsstrahlung as given by Equation [1] For our simu- 
lated electron f(v, x, t) - F(E, x, t)/mf. distribution this will be 



I(e, X, f) - 



4nR2 



fjv, x,f) 
n(x) Q{£, E) 



dEdx 



(6) 



where A is the area of the emitting plasma. We take this to be 
FWHM^ = 8(ln2)t/2, and for the chosen d the area matches 
the characteristic sizes of flaring magnetic loops observed with 
RHESSI, (e.g. Emslie et al. 2003; Kontar et al. 2010). 

The time averaged spectrum from a particular spatial region 
of our simulation (between x\ and xi) is 



E x=x\ t-0 



n(x) 



f(v,x,t) 



Q(e,E) 



dEdxdt 



(7) 



This allows us to calculate the X-ray emission either from the to- 
tal source, coronal source (above 30Mm) and the footpoint (be- 
low 3Mm) as shown in Figure [1] 
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Another quantity of interest is the mean electron flux spec- 
trum (nV FjE)) as it is deducible from HXR spectrum (e.g. 
iBrown et al.. .2006). From our simulations this is calculated as 



{nVF{E)} 



n{x) 



f(v,x,t) 



dxdt 



(8) 



This is model independent value which can be inferred directly 
from HXR observations. 



3. Simulation Results 

3.1. Continuous Injection 

The first situation to consider is that of the steady-state 
continuous injection of the electron beam. This allows the 
most direct co mparison to previous work on purely collisional 
( Brown 1971h and co l hsion al plus wave-particle interactions 
( Hamilton & PetrosianI Il987h transport, though in both cases 
fliey consider a spatially independent solution. In our simulation 
we continuously inject a beam (using 5b = 5 and n\, = lO^cm"-' 
in Equation |5]l at the coronal source until a steady state is 
reached, after about 1.5 seconds. So observed HXR variations 
on timescales less than this are likely due to non-stationary in- 
jection. The resulting steady-state mean electron flux {nVF{E)) 
and X-ray spectrum /(e) are shown in Figure |2] 

For the coronal source emission the spectra are virtually 
identical, which should be expected since the transport processes 
should have little opportunity to act on the distribution at this 
source location. Fitting a power-law over 20 - 100 keV we ob- 
tain the expected 6cs ~ 5 = 5b in the mean electron spectrum 
and yes ~ 6 = 5b + 1 in the X-ray spectrum for the coronal 
source. We choose this energy range to fit over as it is approx- 
imately over this range that the power-law fit is performed in 
HXR imaging spectroscopy. At lower energies the emission is 
predominantly thermal and at higher, background noise. For the 
footpoint there is a slight difference between the electron beam 
(collisional transport only) and beam and waves (collisions and 
wave-particle interactions) simulation below about 30 keV but 
this is very minor. This confirms the previous results that sug- 
gested that the addition of wave-particle i nteractions would have 
little n oticeable affect on the final spectra (i Hamilton & PetrosianI 
Il987h . We also achieve a spectral index of ypp ~ 4 = 5b - 1 , re- 
sulting in Ay a; 2, both as predicted. The total source spectra are 
slightly softer than the thick- target prediction but that is because 
they include a faint coronal source as well. 



3.2. Transient Injection 

The second set of simulations have the same set up as the pre- 
vious case 33.11 except that we only have the injected electron 
beam at the start and we run the simulation until the electrons 
have completely propagated through the system, losing energy 
and thermalizing to the background distribution and leave the 
numerical grid. This case is therefore referred to as a transient 
initial injection. The mean electron flux and X-ray spectra for 
these simulations, integrated over the 1 second of simulation 
time it takes the electrons to completely leave, are shown in 
Figure [3] For the collisional transport case (beam only) there 
are only minor differences with the previous steady-state sim- 
ulations. The magnitude of the simulation spectra are lower in 
this case but that is because previously the electrons were being 
continuously pumped into the system where as here it was only 
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Fig. 6. The X-ray spectrum of an transiently injected electron 
beam simulation over 1 second, for the footpoint (green), coro- 
nal (black) and total (purple) source, with initial spectral index 
of 5b =5 (left) and beam density Mb - 10*^ cm"-'. This is in 
comparison to the case of rib - cmT^ shown in Figure [3] 



done initially. The other difference is that the simulation spec- 
tra are all steeper (softer) than predicted by both the thin- and 
thick-target models. These however are steady-state models and 
we have simulated the temporal evolution of the electron beam 
so would expect some differences. Although we still obtain a dif- 
ference in spectral index of Ay ^ 2 for the collisional transport 
case. 

With the addition of wave-particle interactions resulting in 
Langmuir wave growth we get a dramatically different spectrum 
for the footpoints (and hence total source as well), considerably 
flatter than that found in the purely collisional case. As the coro- 
nal source is again similar in both cases we obtain a difference in 
the spectral index bigger than the collisional case, the addition of 
wave-particle interactions producing Ay a; 3. In this case the X- 
ray spectrum does start to deviate from a power-law in the range 
over which we are fitting, in particular starting to flatten at lower 
energies. This appearance of a break in the footpoint spectrum is 
consistent with those observed in HXR spectroscopy but as we 
are comparing our simulation to HXR imaging spectroscopy in 
which the energy resolution and noise make it unsuitable to fit 
anything other than a single power-law (Battaglia & Benz 200^, 
we will continue to use a single power-law. Additionally the ap- 
pearance of a break is exacerbated by the sharp drop off in the 
spectra at the highest energies, beyond our fitting range. This is 
due to the upper velocity cut-off in the initial electron beam and 
the simulation grid and we will return to discuss this in detail in 

EH] 

The difference between the two cases results from consid- 
erably more wave growth at higher velocities in the transient 
evolving case. This difference is shown in Figure |4] where we 
have plotted the profile of the electron and the wave spectral 
energy distribution at the height of 10 Mm in the simulation 
for both the steady-state case, and for various times during the 
evolving case. From the quasilinear term in Equation we see 
that we get wave growth from a positive gradient in the electron 
distribution df/dv > 1 as well as from the spontaneous emission 
where there is a concentration of electrons, Sf. In the evolving 
case, the electrons with the highest velocities move away from 
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Fig. 4. Electron /(v) and spectral wave energy W{v) distributions for height of lOMm above the photosphere. The dotted line 
indicates the distributions for the continuous injection case once a steady-state equilibrium has been reached (shown in Figure |2]l 
and the solid lines show the time evolution (different colours for different times as indicated) for the transient initially injected 
electron beam (shown in Figure |3)- 
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Fig. 5. The X-ray spectrum of an transiently injected electron beam simulation over 1 second, for the footpoint (green), coronal 
(black) and total (purple) source, with initial spectral index of 5b = 3 (left) and cJb = 7 (right). This is in comparison to the case of 
(5b = 5 shown in Figure |3] 



the bulk of the distribution creating a positive gradient on the 
leading edge resulting in strong wave growth. At lower veloci- 
ties the emission is dominated by the spontaneous emission from 
the bulk of the electron distribution but this has little effect on 
the X-ray spectrum as it is the higher velocity electrons that can 
greatly change the X-ray emission. If this is allowed to reach 
a steady-state with a continuously supply of electrons behind 
those that move away, the wave growth is highly diminished at 
high velocities, since the strong wave growth from the leading 
edge of the electron distribution is washed out by the following 
electrons. This higher rate of wave growth in the transiently in- 
jected case for the high energy electrons provides an additional 
mechanism for them to lose energy to the background plasma 
and hence flattens both the electron and X-ray spectra. 



3.2.1. Initial Spectral Index 

In Figure |5] we show additional simulations of the transient ini- 
tially injected electron beam with different initial spectral indices 
5b = 3 (left) and i5b = 7 (right), instead of 5b = 5 shown in 
Figure|3] For the harder (flatter) spectrum (6^ = 3) there is a ma- 
jor difference between the beam only and the beam and waves 
simulations in both the spectral indices and magnitude of emis- 
sion. This is again due to there being more high energy electrons 
which have an additional energy loss mechanism due to the wave 
growth. 

In contrast, with a steeper (harder) initial spectrum (5b = 7) 
there is a smaller difference between the two transport cases, as 
there are fewer high energy electrons, so there is considerably 
lower rate of wave growth and consequently collisional effects 
dominate throughout. As a result the increase in the difference 
between in the spectral indices between the coronal and foot- 
point sources for the beam only and beam and waves simulations 
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Fig. 7. Difference in the spectral index of the coronal source and footpoint for the mean electron spectrum Ad (left) and X-ray 
spectrum Ay (right) as a function on the spectral index of the injected electron beam 6b and beam density (blue 10^ cm"^, red 
10** cm""*). Shown are the results for the beam only propagation (subject only to Coulomb collisions) and the beams and waves 
simulation. The grey points indicate the results when a higher maximum velocity is used, discussed in ^3.3.11 



is almost half that with the flatter spectrum than with the steeper 
spectrum. 

3.2.2. Initial Electron Beam Density 

In Figure|6]we have increased the initial beam density by an or- 
der of magnitude compared to Figure |3] to = lO^cm"^. For the 
beam only case, the increase in the number of accelerated elec- 
trons has only increased the magnitude of the X-ray emission, 
with the spectral shape, and hence difference in spectral indices 
between the coronal and foopoint sources, identical to the lower 
beam density simulation (right Figure |3]l. 

With the inclusion of wave-particle interactions in the beam 
and waves case we not only see a change in the magnitude of 
emission but also in the spectral shape, with the footpoint source 
flatter than with a lower beam density. This is again due to there 
being more electrons a higher energies, producing higher lev- 
els of wave growth which quickly flatten the emission. This ad- 
ditional energy loss mechanism also explains the difference in 
magnitude of the emission between the beam only and beam and 
waves. It is important to note that the wave-particle interaction 
is a highly non-linear process and so can dramatically change 
the HXR spectra with only a slight change of the electron beam 
parameters. 

3.3. Difference in Coronal & Footpoint Spectral Indices 

In Figure |7] we show the difference in the spectral index of the 
coronal and footpoint source in the mean electron flux (A6) and 
X-ray (Ay) spectra as a function of the initial electron distribu- 
tion spectral index 6b. Considered are both the beam only (col- 
lisional transport, dashed lines) and beam and waves (collisions 
and wave-particle interactions, solid lines) for an initial beam 
density of «b = 10^ cm"-' and nb = lO'^cm"-'. For the purely 
coUisional case we see that there is no difference as the beam 
density nb increases and only a small change with increasing the 
initial spectral index 6b. This change is due to the upper cutoff 
in electron velocity (discussed further in ^3.3.11 ) and the energy 
range we choose to fit over. Again we are choosing this fixed 
energy range as this is comparable to the range over which a sin- 



gle power-law is typically fitted in HXR imaging spectroscopy. 
Note that we would not expect the spectral difference to be ex- 
actly Ay = 2 as this is predicted from the thin and thick-target 
approximations which is closer to our steady-state simulations 
the transiently injected ones presented here. 

There is considerable change however with the addition of 
wave-particle interactions in the beam and waves cases. The dif- 
ference in the HXR spectral indices is consistently larger than 
when wave-particle interactions are included to the purely col- 
lisional case. For an increased beam density there is a bigger 
difference in the HXR spectral indices for 6b > 4, typically an in- 
crease of around 0.7. As the initial spectral index decreases, the 
difference between the different beam densities disappears (at 
6b = 4) and then swaps over for the flattest spectrum. The reason 
for this is that these are the simulations with the largest propor- 
tion of high energy electrons in them and hence the wave growth 
is having a major effect. However at these energies we starting to 
approach the maximum velocity of the simulation grid, so the re- 
sults are dominated by edge effects, which we discuss in 33.3.11 

3.3.1. Higher Maximum Velocity 

To investigate the role of the upper cutoff to the initial elec- 
tron distribution we reran the most severely affected case with 
the most electrons initially at high energies {6b - 3 and nb = 
10**cm"-') with a higher maximum velocity of vq - 200vt (£0 = 
1730keV), compared to vq = 75vt (£0 = 243keV) used pre- 
viously. The resulting mean electron flux and X-ray spectra for 
these simulations in comparison to the previous ones are shown 
in Figure |8] In the mean electron spectra there are no changes 
in the 20 to 100 keV energy range we have been fitting with a 
power law for the coronal source or in the collisional case. In the 
footpoint source spectrum when wave-particle interactions are 
included there is a slight dip about 50 keV in the spectrum which 
extends to higher energies. The major changes appear in the X- 
ray spectra as these are integrated over all electron energies, as 
can be seen in Equation [T] The footpoint source in particular is 
badly affected, with the higher cutoff producing a spectral index 
that is 0.87 flatter. 

Returning to Figure Q we can see the effect of this higher 
maximum velocity in the difference in the coronal and footpoint 
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Fig. 8. Mean electron flux spectrum (left) and X-ray spectrum (right) of the transiently injected electron beam simulation over 
1 second, for the footpoint, coronal and total source. Shown are the simulations using a maximum velocity of vo = 75vt (Eq - 
243keV) (black coronal source, green footpoint) and vq = 200vt (Eq - 1730keV) (grey lines). 




Fig. 9. .(Right) The time averaged beam energy U(x) for energy bands 7-20 keV, 20 - 50 keV and 50 - 100 keV (blue, orange 
and green respectively) for simulation with collisional transport (dashed lines. Beam Only) and additional wave-particle interactions 
(solid lines. Beam & Waves). (Left) The beam energy loss rate, normalised by the beam energy, for the same energy bands and 
simulations shown in the left panel. 



spectral indices (grey points). In the mean electron flux distri- 
bution there is only a minor change but there is a substantial 
change in the X-ray spectral indices. With the collisional case 
the difference in spectral indices returns closer to Ay x 2. With 
the wave-particle interactions we get a bigger change in the X- 
ray spectral indices 2 < Ay < 4, similar to X-ray observations. 



Given the issue of the maximum velocity cutoff the natural 
question then arises over the non-relativistic nature of our simu- 
lations. We do not feel that we need to consider a relativistic ver- 
sion as we are tryin g to make a di rect comparison to the collision 
thick-target model (iBrownlfTOTlh which itself is non-relativistic. 
Also the effect of relativism on electrons around 100 keV is min- 
imum (more importantly it is much less than the effect of upper 
energy-cutoff) and wave-particle interactions are weaker above 
100 keV. 



3.4. Beam Energy 

The energy content of the beam is calculated as 

U(x,t)^^J^'v^-f(v,x,t)dv, (9) 

so variations in the beam energy at a particular time or posi- 
tion is due either to there being a change of electron velocity 
or number of electrons in that velocity/energy range. This can 
either occur from the particles gaining or losing energy to the 
background plasma or the beam moving through the particular 
region. In Figure |9] we have plotted this beam energy, time aver- 
aged over the whole 1 second of simulation time and normalized 
by the maximum energy, for various energy bands and both pro- 
cesses. Here we see that the simulations with wave-particle in- 
teractions quickly lose energy in the corona. This again matches 
what we have found with the flattening of the mean electron and 
X-ray spectra, for example Figure |3] where the waves provide 
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Fig. 10. . (Left) HXR spectra of the footpoint region in the simulation with 60 = 5, «b = lO^cm"^ for various starting heights 
xo - 10, 40, 60 Mm (green, orange and blue respectively) for both purely collisional and beam and wave interactions. (Right) The 
fitted spectral index over 20 to 100 keV of the footpoint region as a function of electron beam starting height xq for purely collisional 
transport (beam only) and additional wave-particle interactions. 



the electron beam an additional energy loss mechanism to the 
background plasma , heating the leg of the loop instead of the 
footpoints. In the purely collisional case, the energy content of 
the beam between 20-50 keV and 50-100 keV is 55% and 90% of 
its initial coronal value when it reaches the chromosphere, only 
45% and 10% lost to heating the leg. With the additional wave- 
particle interactions the beam in these energy ranges losses the 
majority of its energy in the corona, with 90% and 80% of its ini- 
tial energy being lost to heat the loop at heights above 20 Mm. 

We can use this time averaged beam energy to calculate the 
energy loss rate as a function of spatial position dU/dx, shown 
in the righthand side of Figure|9] For the Coulomb collision only 
simulations {Beam Only) we see that the energy loss rate closely 
matches the profile of the background plasma density: it is con- 
stant in the corona and sharply rises through the transition region 
and chromosphere, i.e Figure[T] For the cases where we have in- 
cluded wave-particle interactions we see a similar structure to 
that of purely collisional transport only high in the corona and in 
the chromosphere. In these regions the waves have not been able 
to grow and the electron transport is dominated by collisional 
losses respectively. In the rest of the corona the Beam &■ Waves 
case deviates from the purely collisional simulations, the energy 
loss rate being far higher for the majority of the beam's transport 
in the corona. 

3.5. Starting Height: Difference in Footpoint Spectral Indices 

Another parameter of the initial electron beam to vary is the 
starting height xq. This is of interest as it might be able to ex- 
plain the observed difference in spectral index between double 
footpoint sources. In a symmetrical injection there would be the 
same starting height (or path length) from acceleration/injection 
site down to the footpoint. In an asymmetrical loop one foot- 
point would be closer to the site of injection and the other fur- 
ther away. This can be approximated in our simulations by using 
different starting heights given that the coronal background den- 
sity profile is the same. The results for various starting heights 
are shown in Figure [10] For all heights the result for purely col- 
lisional transport is very similar but varies when wave-particle 
interactions are considered. For the lowest starting height there 
is little time for wave growth to develop to any significant level 



and so the resulting electron and HXR spectrum is similar to 
the purely collisional case. So for the initial beam parameters 
shown we can achieve a difference in footpoint spectral index 
of around Aypp ~ 1 from just changing the starting height of 
the injection. This would have the largest effect for a loop with 
one leg path length from the acceleration site to footpoint of 
< 20Mm and the other > 20Mm: either a short 30Mm with 
a small amount of asymmetry or a highly asymmetric longer 
loop. Typicall y for microflares the SXR thermal loop is typically 
30Mm long ( Hannah et al. I I2OOI though this would be shorter 
than the combined distance from the acceleration site to both 
footpoints. The time of flight distance (assumed acceleration site 
to HXR footpoint) was f ound to range between 10 and 50 Mm 
dAschwanden et al.|[T996l) . 



4. Discussion & Conclusions 

The simulations presented here show the importance of in- 
cluding wave-particle interactions and to simulate the time and 
spatial evolution of flare emission. The growth of Langmuir 
waves flattens the flare spectra producing a larger spectral dif- 
ference between the coronal and footpoint sources compared to 
purely collisional transport. This might be able to explain the 
Ay > 2 observed in some flares with RHESSI. Also, when 
considering asymmetrical loop injection, wave-particle interac- 
tions maybe able to explain the differing spectral indices be- 
tween pairs of simultaneous footpoints without modifying the 
background plasma density profile. Langmuir waves is the dom- 
inate process in the corona and provides and additional mech- 
anism for the beam to quickly lose energy to the background 
plasma. This however does result in fewer high energy elec- 
trons and fainter HXR spectra than those from purely collisional 
transport. Note that they quickly lose their energy to a few keV 
and then take longer to reach closer to the background ther- 
mal level. The wave-particle interactions therefore help produce 
spectral indices closer to observations but require a larger num- 
ber of injected electrons in comparison with collisional trans- 
port. Ways to alleviate this problem could include additional 
chromospheric energisation (Brown et al. 2009) or wave refrac- 
tion from inhomogeneities in the background plasma (iKontail 
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I2001al: iReid & Kontdl20T(l . Both are currently under investi- 
gation and subject to future publications. 

This work is a step towards a more complete treatment of 
electron transport in solar flares with many other processes still 
to be considered. For instance we do not consider the anisotropy 
of the electron beam or the convergence of the magnetic field 
and simulate only one spatial dimension. The Langmuir waves 
generated in our simulations are the seeds for radio emission 
which are sometimes seen in reverse type 111 bursts. Therefore, 
a natural next step for this work would be to calculate the radio 
emission from the propagating electron beam. This is a compli- 
cated task given the highly non-linear processes involved in the 
emission. By using both radio and highly sensitive HXR obser- 
vations in comparison to simulations such as those presented in 
this paper we would be able to greatly constrain the nature of 
particle acceleration and transport in solar flares. 

The HXR imaging spectroscopy of RHESSl has shown the 
importance of such observations for solar flare physics. However 
to really pin down the acceleration and transport processes 
would require having solar HXR imaging spectroscopy instru- 
ments providing orders of magnitude better sensitivity and dy- 
namic range over RHESSl, so that the coronal source emis- 
sion is routinely observed with the bright footpoint source. Such 
advances in sensitivity and dynamic range are a chievable with 
HXR focusing optics such as those on FOXSI ( iKrucker et alJ 
[2009) and NuSTAR (Harrison et al. 2010). 
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